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ABSTRACT 

Aims. Accurately measuring the atmospheric coherence time is still an important problem despite a variety of applicable methods. 
The Multi-aperture scintillation sensor (MASS) designed for the vertical profiling of optical turbulence, also provides a measurements 
of coherence time, but its results were found to be biased. Hence there is a need for a more robust method to determine To. 
Methods. The effect of smoothing the stellar scintillation by a finite exposure of the detector is considered. The short exposure regime 
is described and its limits are defined. The re-analysis of previous measurements with the MASS is performed in order to test the 
applicability of this approach in real data processing. It is shown that most of the actual measurements satisfy the criteria of short 
exposures. 

Results. The expressions for the mean wind speeds V 2 in the free atmosphere from the measurement of the scintillation indices 
are derived for this regime. These values provide an estimate of the atmospheric coherence time To without the need of empirical 
calibration. The verification of the method based on real measurements of the resulting r are in good agreement with independent 
methods. 

Key words. Site testing - Instrumentation: adaptive optics - Atmospheric effects 



1. Introduction 

In recent decades the adequacy of optical turbulence (OT) mea- 
surements in the Earth's atmosphere above potential astronomi- 
cal sites and operating observatories has become more and more 
important. Developments in astronomy demand more in the way 
we determine OT. One such requirement is the ability to use 
cheap tools to monitor the OT on long time scales and automat- 
ically. 

For such methods using a modest feeding optics, the analy- 
sis of the spatial wavefront distortion is replaced by the anal- 
ysis of temporal variations of the wavefront inside a lim- 
ited are a. The widespread Differe ntial Image Motion Monitor 
DIMM dSarazin & Roddien [1990) which measures the turbu- 
lence integrated over the whole atmos phere, and the Multi- 
Aperture Scintillatio n Sensor MASS dKornilov et all [2003; 
Toko vinin et al.l20Tj3h which determines the altitude distribution 
of the OT belong to such methods. 

The transition from a spatial to a spatial-temporal desc ription 
of the OT is based on the hypothesis of frozen turbulence (Tavlo 
1938) a nd wa s analyzed by many authors (see e.g. lRoddierll981 
Martini [T987I) . A necessary component to this analysis is the 



knowledge of wind velocities in the atmosphere. 

The MASS instrument was developed to measure the verti- 
cal distribution of the OT. Its data contain the necessary infor- 
mation about the temporal properties of stellar scintillation to 
estimate the atmospheric coherence time tq with using a method 
involving the differ ential exp osure scintillation index DESI as 
described in Toko vinin! (120021) . 

After the actual measurements of OT with the MASS, it 
became clear that to values derived in this way were signifi- 
cantly underestimated, and therefore re quired additional calibra - 
tions and corrections during analysis dTravouillon et al.l 12009). 
Empirical recipes of correction do not give confidence results 



because this systematic u nderestimation may be due to several 
factors (lTokovininll201 ll) . On the other hand, a large and per- 
manently increasing volume of collected data requires an uni- 
form treatment for objective comparison of OT characteristics 
between various sites. 

In this paper we analyse the effect of temporal averaging 
which affects on scintillation indices measured in the MASS in- 
strument and their role in the evaluation of to. In addition the 
modification of the DESI method which does not depend on an 
empirical calibration is proposed. The presentation of the modi- 
fied method is preceded by theoretical description of scintillation 
in the regime of short exposures. 

The regime of short exposures allows us to establish a sim- 
ple relationship between the measured scintillation statistics and 
the characteristic wind speed entering the definition of to. In the 
last section, the meth od is tested with MAS S/DIMM data taken 
at Mt. Shadzhatmaz dKornilov et alj 1201 Oh and Mt. Maidanak 
dKornilov et al.ll2009h . 

We found that the main cause of the mentioned underestima- 
tions of To was the wrong interpretation of the formula for DESI 
and these discrepancies virtually disappear after the correction 
in the data processing software. 

The modified method does not require any empirical calibra- 
tion. It is theoretically clearer, and gives an estimate of the mean 
wind in the free atmosphere as well as a more accurate atmo- 
spheric coherence time. 



2. Temporal averaging of the scintillation 

The theory of weak perturbations implies that the scintillation 
index s 2 - the variance of relative fluctuations of light intensity, 
is described by the sum of the scintillation indices produced by 
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individual turbulent layers: 



o 



C 2 n {h) W(h)6h, 



(1) 



where W(h) is the weighting function (WF) which depends on 
the size and shape of the receiving aperture and does not depend 
on the altitude distribution of the structural coefficient of the re- 
fractive index C 2 (h). WF represents a power of the scintillation 
generated by a layer of unit intensity located at a height h. 

The MASS method involves simultaneous measurements of 
the scintillation indices in 4 concentric apertures of different 
diameters (hereafter: A, B, C and D), leading to 10 indepen- 
dent scintillation indices. The vertical OT profile is restored 
from the measured indices and the theoretic ally calculated WFs 
dTokovinin et alJl2003b iKornilov et alJl2007h . The calculation of 
a set of functions W(h) is based on the assumption that light in- 
tensity measurement has a zero exposure, i.e. the averaging fac- 
tor is related only to the receiving aperture. 

The scintillation index s 2 measured with finite exposure time 
is determined by the integrated effect of all turbulent layers along 
the line of sight as in the case of zero exposure ([1} but with 
another WFs. 

The expression for the new WFs W(w, t, h) depends not only 
on the layer altitude h but also on the wind speed w = w(h) 
and averaging time r dTokovininl2002tlKornilovl201 lb . It differs 
from the expression for W(h) with an additional multiplicand 
A s {w, t, /) in integrand: 



W'(w, t, h) = 9.62A- 2 f- 8/3 sin V 2 f) A(f) A s (wrf) df. 
Jo 

(2) 

Here, / is the modulus of the spatial frequency, A(f) is spa- 
tial aperture filter (axisymmetric for MASS apertures), rf is the 
Fresnel radius r 2 F = Ah. 

Multiplicand A s (wt/) = A s (0 is the spectral filter averaging 
with the wind shear wt and can be expressed through t he Bessel 
funct ions Jq and J\ and Struve functions Ho and Hi dKornilovl 
1201 lh : 

AM) = 2J (2nO-^^-nM2nOH 1 (27T{)+nJ 1 (27:f)H (2n{). 

(3) 

It has a simple asymptotic behavior for small values of argument 
£ = wr/: A s (0 « 1 - 7r 2 £ 2 /6, as it follows from its series ex- 
pansion in the neighborhood of 0. The approximation provides 
an accuracy better than 0.02 until < 1. When £ — » oo, the 
function A s (£) m \jn^, and starting from£ m 1, the relative dif- 
ference is less than 0.04. 

These asymptotes correspond to the two extreme cases: the 
short and the long exposur e regimes. Regime of short exposures 
(SE) considered earlier by Tokovinin (2002) will be the subject 
of this pa per. Regime of long exposures (LE) was studied re- 
cently bv lKorniiovl d201 lh in the application of potential impact 
of stellar scintillation on the accuracy of photometric measure- 
ments. 

The main feature of these regimes is that the weighting func- 
tion W(w, t, li) can be represented as the product of a function 
depending on wind shear, and the function which is independent 
on winds. It allows to separate the wind effect and the geometry 
of light propagation. 



3. Short exposures 

In the expression (f2]i the integrand is significantly different from 
in the region of intersection of the aperture (/ < l/D) and 



Fresnel (/ < 1 /rf) filters. Outside this domain, i.e. when / > 
minj 1 /D, 1 /rf} the integrand tends rapidly to 0. If the wind shear 
tw <K max{D, rf } then frw < 1 and hence a quadratic approx- 
imation 1 - 7t 2 (wt/) 2 /6 of the spatial filter A s (w, r, /) can be 
applied. 

Let's replace A s for a certain aperture with the quadratic ap- 
proximation in the expression (f2} and take w and r outside the 
integral over frequency. Then 



s 2 = j C 2 „(h)W'(w,T,h)dh = 

=! c - 



Ah) W(h)dh- 



2 (h) w(h) 2 U(h) Ah, (4) 



In this difference the first integral is the scintillation power s^ 
for zero exposure. The second integral (we denote it as *VF) is 
the atmospheric second moment of wind additionally weighted 
with U(h). Weighting function U(h) is obtained by multiplying 
the initial spatial spectrum by 7r 2 / 2 and subsequently integrating 
over /. Since the effect of high-frequency spectral components 
increases after such multiplication the functions U(h) become 
essentially different from the usually used W(h) and have dimen- 
sions of irT 7 / 3 . The set of U(h) functions is shown in Fig.Q] 
Additional features of these WFs are as follows: 

- The spectral band width effect is more important than before 
owing to the high-frequency spectral components which are 
more intense. The difference of WFs for white and red stars 
reaches 10 % in aperture A. 

- For infinitely small aperture an asymptote does not exist 
since at D — the integral diverges. However, there is an 
envelope of ~ h~ l t 6 to which the curves converge as D — > 
(see Fig. |T}. 

- The asymptote for the aperture D » rf is also very interest- 
ing: U(h) * 17.22 ,T 2/3 D 3 h 4/3 . If compared to the asymp- 
totic behaviour of the normal scintillation index, it is chro- 
matic and has a stronger dependence on D. 

Note, that the WFs shown in the figure by the dashed lines 
are calculated for cross-indices, i.e. the values obtained directly 
from the MASS instrument measurements rather than differen- 
tial indices. Differential indices were always calculated from the 
normal indices and covariances in the form of linear combina- 
tions, but from the viewpoint of the OT restoration it is a super- 
fluous intermediate step which we have dismissed. 

The SE regime and the LE regime, are self-reproducible in 
the sense that if they occur in each turbulent layer they will be 
observed for the whole atmosphere. To determine the limits of 
SE regime, one can use the condition of applicability of the 
quadratic approximation of the wind shear filter A s (%), namely 
that at its border the function is equal to 1 - 7r 2 £ 2 /6 = 5/6. 

Since for small wt this function varies slowly with frequency 
/ within the main spectral peak (its maximum is located at fre- 
quency < 0.7 /rf) of the integrand, then the value of s 2 will 
also slightly vary with wt. Consequently we can assume that 
the quadratic approximation of the form (|4} for a given aperture 
and atmospheric conditions is applicable while t 2 "VF < s^ or 



2 5 9 
«r > ^ 



(5) 



Of course this inequality defines the SE regime only within 
some level of accuracy. One should take into account the error 
of approximation of the function A s (£) (« 0.02) and the contri- 
bution of the secondary peaks of the scintillation spectrum. A 
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Fig. 1. Weighting functions U(h) for a set of MASS/DIMM aper- 
tures computed for a typical spectral response of its detectors and 
the light source of spectral class AO V. Solid curves - WFs for 
normal indices, dashed - for cross-indices. Thin lines denote the 
functions for a number of smaller apertures, their diameters are 
marked in centimeters 




4 8 16 
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Fig. 2. Indices s 2 normalized to as function of an exposure t. 
Bold curve corresponds to the median of 712 = 0.938, dashed 
line - to the limit of the SE regime 712 = 0.870. Another curves 
are marked by appropriate y\%. The upper horizontal line repre- 
sents the lower boundary of SE regime, the lower horizontal line 
- the upper limit of the LE regime 

Table 1. Characteristics of 721 and 741 distributions for measure- 
ments at Mt. Shadzhatmaz in May 2009 - June 2010 period 



rough estimate of the accuracy of the numerical factor in (|5) is 
« 0.05. 

Similarly it is possible to write the limit of the LE regime in 
the form of s 2 < -s\. Both boundaries are indicated in Fig. [2] 
with dashed lines. The Figure demonstrates examples of typical 
behavior of the measured indices s 2 with exposures ranging from 
1 to 64 ms in different wind conditions. Note that the left and 
right branch of each curve can be generated by different turbulent 
layers because in SE regime the contribution of an individual 
layer is ~ C 2 w 2 while in the LE regime ~ C 2 /w. 



Parmeter 


A 


BCD 


AB 






Distribution of 721 




Median 


0.94 


0.95 0.96 0.97 


0.95 


Out of SE regime, % 


18.0 


12.5 6.3 3.3 


13.8 






Distribution of 741 




Median 


0.77 


0.80 0.84 0.87 


0.79 


Out of SE regime, % 


67.0 


59.6 49.8 39.6 


61.3 



4. Verification of applicability of SE regime in the 
MASS measurements 

Since the value is calculated and not measured the restriction 
(f5]l is not directly applicable. However instead of one can use 
an additional index measured with a shorter exposure r' < t. 
Simple algebraic manipulations lead to the condition 



It follows that for exposures t — 2 ms and t' = 1 ms s\ should 
be more than 0.870 s\, for 3 and 1 ms s\ > 0.850 s 2 and for 4 
and 1 ms s 2 A > 0.842 s\. Condition © can not be used when 
exposure r' itself is too long for the SE regime. 

The analysis of actual MASS/DIMM data obtained at 
Mt. Shatdzhatmaz in the period May 2009 to June 2010 was 
performed with the usual exposure r = 1 ms. The special ver- 
sion of MASS software recorded additional data in output files 



to calculate later the indices with exposures 2, 4, 8, 16, 32 and 
64 ms (see Fig. |2j. 

Cumulative distributions of ratios 721 = and 741 = 

s\l s\ averaged over a minute are shown in Fig. [3] As it might 
be expected from the form of the WFs (Fig. Q3 smaller ratios are 
observed for smaller apertures, more sensitive to OT motion, and 
the distribution for AB cross-index is very close to the curve for 
the aperture B. Numerical characteristics of the distributions are 
given in Table Q] and show that the 2-ms measurement is almost 
always within the SE regime. Only 18% of the measurements in 
A aperture fail to meet the 721 > 0.870 condition. 

Characteristics of the distributions can vary greatly from sea- 
son to season. E.g. during February-March 2010 period, which 
was characterized by strong winds, the median 721 decreased to 
0.90, while the fraction of measurements that do not satisfy the 
criterion of SE increased to 38%. On the other hand, all measure- 
ments fell into SE regime in October 2009 which was notable for 
its stable calm weather. 
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Table 2. Characteristics of 721 and 731 distributions for measure- 
ments at Mt. Maidanak in August 2005 - December 2007 period 



Parmeter 



B 



C 



D 



Median 

Out of SE regime, % 
Median 

Out of SE regime, % 



Distribution of 721 
0.92 0.93 0.94 0.96 
19.4 13.3 5.5 0.3 

Distribution of 731 
0.82 0.84 0.86 0.91 
64.1 57.6 41.5 12.8 




0.995 0.99 



Fig. 3. Cumulative distributions of 721 (black curves) and 741 
(grey curves) for the apertures A, B, C, D and AB cross-index 
based on measurements at Mt. Shatdzhatmaz. Vertical lines cut 
off points which do not satisfy the SE regime 



For measurements with 4 ms exposure (value 741) the situa- 
tion is radically different. In this case, the vast number of mea- 
surements does not satisfy the SE condition and consequently 
the quadratic approximation can be used only in rare cases and 
with caution. 

Similar analysis was performed for measurements at 
Mt. Maidanak. The standard MASS data were used, from which 
indices s\ = (s\ + pi)/2 and s\ = (3s 2 + 4pi + 2p2)/9 were 
calculated to construct cumulative distribution of 721 and 731 = 
s 2 /s 2 . General properties of these distributions do not differ from 
that of Mt. Shatdzhatmaz. Some differences are due to signif- 
icantly larger aperture s C and D of the first generation device 
dKornilov et alj [2003). Characteristics of the distributions are 
given in Table [2] The fraction of measurements dropping out of 
SE regime in the case of A aperture is 19.4% for exposure 2-ms 
and 64% for 3-ms exposure. 

In the situation where 721 approaches unity, a very high pre- 
cision of the ratio is required as the curves for 721 = 0.995 and 
721 = 0.99 differ very much. For these curves even 8 ms mea- 
surements are in the SE regime. 



5. Reduction to zero exposure 

Typical exposures of 1 ms are taken with the MASS. This leads 
to a wind shear of the order of 3 cm assuming wind speed in 
the tropopause ~ 30 ms -1 . The value is comparable with the 
size of the device apertures and with typical Fresnel radius so 
1 ms exposure can not be considered as infinitely small one. The 
procedure used to correct indices to zero exposure is provided 
during the MASS data processing. Th e algorithm is based on 
numerical simulations dTokovin in 2002). 

However for the measurements in SE regime, the required 
correction can be calculated with the help of a direct method 
using the two indices s\ and s\ obtained with different exposures. 
Although the ratio of exposures may be arbitrary, it is convenient 
to consider the case of a single r and a double 2t exposures. 
Then the expression (|4]l becomes: 



- < vk 



After solving this system for s2 we obtain 



4r 2 



-<VX 



(7) 



(8) 



One can express s\ in terms of s. and covariance p\ of adjacent 



counts as 2s\ 



s, + p\. The formula can be re- written as: 



¥ 1 



6 Pl - 



(9) 



The resulti ng correction i s som ewhat smaller than the one 
adopted in iKornilov etafl d2007l) : s 2 = 1.25s 2 - 0.25pi. Of 
course, formulas (|8} and (|9) are exact only if the condition (fSJ) of 
SE regime is satisfied for s\. 



The statistical error of 



Pi 



very 



close to that of s, 



Consequently the correction of the scintillation index increases 
its standard error by a factor of 1.18 which is not significant in 
practice. 

6. Atmospheric second moment of wind 

We can estimate the integral T^ 7 from the system of equations 
(0. Solving the system for the indices measured in the /'-aperture 
(j = I,..., 10) with exposures T\ and T2 relative to this un- 
known, we obtain: 



2 (h) w(h) 2 Uj(h) Ah : 



(10) 



where Aj denotes the measured quantity. The integral is of little 
interest by itself because it includes an additional factor UjQi) 
which distorts the contribution of different altitudinal layers. 
Real estimation of the second atmospheric moment of the wind 
% = / C 2 (h)w(h) 2 dh can be obtained if we find a linear com- 
bination A u(h) of functions UjQi) which is close to unity. 

Such combination is shown in Fig. [4] Good results can even 
be obtained when only using functions corresponding to normal 
indices. The method and results of the decomposition A u(h) « 1 
are described in detail in Appendix A. 

Expansion coefficients of Au(li) = Yjj c jUj(h) are of order 
of 10" 15 m 7 ^ 3 . Summing equation (fTUl i for the indices with these 
coefficients we obtain: 



(11) 
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Fig. 4. Approximation A u (h) of a constant level with help of the 
set of U(h). Black curves (1) show the approximation with 10 
indices, grey (2) - with 4 indices. For clarity the grey graphs are 
shifted to the level of 0.8. Thick lines correspond to the weight 
of ~ h, thin - to other power of altitude. Curve 3 corresponds to 
the weighting of ~ /i 7 ' 5 . Gray dashed 4 - WF for the aperture A 



It is clear that the boundary layer (h < 0.5 km) wind will not 
be accounted for fully in the integral and the surface wind will be 
completely excluded despite the strong surface turbulence. The 
behavior of Au(h) is similar to the approximating function un- 
der the integral defining the intensity of turbulence in the free 
atmosphere Jf ree . Therefore this quantity can be used for fur- 
ther normalization. The contribution of the surface layer can be 
accounted for later on base of DIMM data. 

Dividing both sides of the expression (fTTT i by Jf ree , we obtain 
the mean square of the wind speed (w 2 ) in the free atmosphere 
and its expression in terms of measurable quantities: 



<w 2 > 



^2 _ ZjCA] 



Jfre, 



Jfre 



(12) 



7. Atmospheric coherence time 



It was shown by iKellerer & Tokovininl ([2007) that the mean 
wind V2 can be a good estimator of V5/3 which is used in the 
definition of the atmospheric coherence time tq. Recalling the 
fact that the instrum ent FADE (Fast Defocusing of stellar image, 
iTokovinin et alj2008l) measures the time constant (interferomet- 
ric coherence time) also using Although authors indicate that 
on average V2 ~ I.IV5/3, we will assume that both values are 
equivalent in the following argumentation. 

Evaluation of V2 (in the free atmosphere) from the measure- 
ments is obtained directly from the formula for (w 2 ): 



V 2 = (w z ) 1,z = 



2xl/2 _ (XjCjAj) 1 ' 2 



jl/2 

free 



(13) 



Substituting this value in the definition of the atmospheric time 
constant tq we obtain the following expression: 



T = 0.314-5- = 0.058A 6/5 J, 
V2 



-3/5 



./ 



1/2 
free 



(14) 



where we have substituted the Fried parameter ro with OT in- 
tensity J lot in the whole atmosphere. For the wavelength A = 
500 nm: 



T = 1.593 ■ 10" 9 /, 



-3/5 



1/2 
free 



(15) 



If there is information about the OT power in the boundary 
layer and surface wind speed Vq, then its contribution can be 
accounted for at the stage of calculating the mean wind V2. In 
order to do this, we correct the mean square of the wind speed 
by adding V 2 in proportion of the surface layer intensity Jql'- 



<w 2 > 



_ (W 2 ) J f ree + V 2 J CL 



or 



V 2 



{ZjCjAj + V 2 JG L ) 



1/2 



1/2 



(16) 



(17) 



Substituting the corrected value in the formula ([TBI we obtain 
an estimation of the atmospheric coherence time for the whole 
atmosphere: 



T = 1.593 ■ 10" 



1/10 



Note that the usual method of adding coherence time t q 



(18) 



5/3 



-5/3 
l GL 



+ T 



-5/3 
free 



is not correct if used to calculate V2 which has 
square metric. 

The relative accuracy of to does not practically depend on 
the accuracy of the J tot , and is twice as good as the accuracy of 
the measured A,. 



8. Verification of the method of evaluation r 

The method presented here was tested with data obtained at 
Mt. Shadzhatmaz and Mt. Maidanak which were analyzed in 
Sect. [4] To construct ^jCjAj we have used the coefficients from 
the Table IA. ll defined for stars of spectral class AO V. This sim- 
plification introduces a systematic error when processing the 
measurements of stars another spectral classes so further results 
should not be regarded as final (less than 5%, see Sect. |3). We 
also computed the coherence time for free atmosphere, i.e. Jf ree 
value was used instead of J, , in formula (TT~5T >. 

For the measurements at Mt. Shadzhatmaz differences Aj 
were computed from indices s\ and s\. The results of the eval- 
uation are shown in Fig. [5] where the cumulative distributions 
of To obtained with different decompositions are presented. The 
method results in 6.69 ms median when using 4 normal indices 
and 6.53 ms using all 10 indices. The difference between these 
distributions is completely explained by the fact that usage of 
all the indices involves additional contribution of moving turbu- 
lence between 0.5 and 1 km. 

In the case of the standard MASS output (Maidanak cam- 
paign) normal indices s 2 , s\ and s 2 were used to study biasing of 
to estimations with different temporal span. At first we investi- 
gated the effect of the DESI formula correction, i.e. replacement 
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5 10 15 20 25 



t q! ms 

Fig. 5. Cumulative distribution of free atmosphere to from mea- 
surements at Mt. Shatdzhatmaz obtained using s 2 - s 2 differences 
with decomposition in 4 (solid line) and 10 (dashed line) indices 




5 10 15 20 25 



t q! ms 

Fig. 6. Cumulative distribution of free atmosphere to from mea- 
surements at Mt. Maidanak obtained using DESI (grey dashed - 
before software correction, grey solid - after) and using s\ - s\ 
and s 2 - s 2 differences (black solid and dashed curves respec- 
tively) 



of the index si with the index s\. As can be seen from Fig.|6]the 
correction increases to estimations by a factor « 1 .4. The cumu- 
lative distributions of coherence time calculated from s 2 - s 2 and 
s 2 - s 2 are presented in the Figure too. One can see that for the 
case of corrected DESI algorithm (5.82 ms) and 2 ms median 
(5.98 ms), the median differs only slightly. 

The to values obtained with si are systematically larger by 
as 0.6 ms than to calculated with s\. This shift is caused by the 
3 ms exposure beyond the SE regime. The larger to gets the more 
3 ms measurements satisfy the criterion of SE and because of 
that, the relative difference between the curves for 2 ms and 3 ms 
is reduced, becoming less than 4% at To as 15 ms. To obtain the 
systematic shift of the 2 ms curve, additional measurements with 
a shorter exposure are required. However taking into account the 
quadratic dependence of indices on exposure, we can estimate 
that its median is not overstated by more than 5%. 



Random errors of to were computed during the averaging of 
V2 over 1 minute as the standard error of the mean. The relative 
error of V2 is approximately constant in entire range of values. At 
very low winds (less than 5 m s _1 ) errors increase. The median 
of relative error is « 0.02, only 1% of the measurements have 
errors greater than 0. 1 . 

9. Discussion and Conclusion 

The analysis of actual MASS measurements described in Sect. [4] 
has shown that for the commonly used 1 ms exposure 80% of 
cases, are within the SE regime and that the correction to zero 
exposure (0 is suitable. Despite some conventionality of the 
regime boundary definition, it can not be significantly weakened. 
On the contrary the study of real data requires some tightening 
of the limit to guarantee that all turbulent layers are within the 
SE regime. Our statistics shows that exposures of 3 ms rarely 
satisfy the criterion of SE. 

Such a control is important because measurements of to 
in large wind shear will be affected by systematic errors. 
Fortunately such cases are not as interesting as opposite ones. 

The reaction of the DESI method on fast turbulence motion 
is not so evident. One ca n derive from the definition of the dif- 
ferential exposure signal (lTokovininl201 ll) and formula dT0T >: 

DESI = ^(3si - 4pi + P2 ) = -(A* 2 - A' 3 ). (19) 

This means that the DESI method estimates the coherence time 
not by using the atmospheric wind moment but its derivative 
with exposure. DESI differences are smaller by an order of 
magnitude than the values themselves and therefore the DESI 
method is noisier. Indeed, for Mt. Maidanak data the median of 
relative random errors is ss 0.03 using the difference s 2 - while 
using DESI the median of relative errors amounts « 0.1. 

The case of calm atmosphere is most interesting from the 
viewpoint of adaptive optics and interferometry. In this situation 
indices differences are small and poorly defined owing to an in- 
crease of relative error. The most dangerous is a systematic error 
due to an incorrect accounting for photon noise. To reduce the 
impact o f such errors, an adaptable choice of exposure may be 
realized ( Tok ovininl20 lib . 

Fortunately, in real situations this effect is not too large. For 
example measurements at Mt. Maidanak show that error Ap = 
0.01 in the parameter p for photon noise correction changes to 
median by 0.05 ms. At To » 15 ms (95% level of cumulative dis- 
tribution) this error results in shift of the curve by 0.7 ms when 
used s 2 - s\ and 0.3 ms when s\ - si. 

Generally speaking, the SE regime is not very suitable to 
determine wind conditions. Much more accurate results can be 
obtained by the determination of the averaging time t for which 
the index s 2 becomes equal to si/2. Such method is similar to 
the ca lculation of correlation peak width, see e.g.. lCaccia et al.l 
(119871) . However, in the middle domain (see Fig. 13 the behav- 
ior of s 2 indicates that similar integral characteristic are closer 
to mean wind speed (w) rather than to the mean square (w 2 ). 
The mean speed may be used for some applications but for to 
evaluation the value (w 2 ) is more suitable. 

The analyzed data were obtained at sites with rare high- 
altitude jet streams. For other observatories, especially closer to 
the equator, the situation could be worse. In such case, it makes 
sense to reduce the exposure to 0.5 - 0.7 ms, which is equivalent 
to increasing the wind range 1.5-2 times. The main objective 
of a correct choice of exposure is to ensure that we are within 
the SE regime for statistically overwhelming fraction of cases. 
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The main advantage of making measurements in the SE 
regime is the simplicity of theoretical description and its prac- 
tical usage. The described method for to evaluation does not 
contain any empirical calibration, it gives more accurate esti- 
mations and additional information about mean wind in the free 
atmosphere. 

All quantities on the right hand side of ( fTBT l are mea- 
sured with MASS/DIMM simultaneously. However the joint 
MASS/DIMM data processin g was realized only i n the recent 
version of the program atmos dKornilov et alJl2010l) and the al- 
gorithm described in Sect. [7] will be implemented in the next re- 
lease of this version. An alternative way is to use a fairly simple 
post-processing of existing MASS output, but in both cases the 
meteo data on the surface wind is additionally needed to take the 
ground layer into account. Free atmosphere to evaluation using 
both methods is implemented starting with atmos-2. 97.3 version. 
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Appendix A: Decomposition in the weighting 
functions 

The decomposition of the atmospheric wind moment 'V 2 is done 
in the same way as an expansion of altitude atmospheric m o- 
ments with sets of measured indices (Tok ovinin et al.ll2003l) in 
the program atmos. In this case we solve the system of linear 
equations Uc = 1, where U is the WFs matrix with dimension 
k x n, c is the vector of coefficients, 1 is the unit vector corre- 
sponding to the altitude grid h. The number of nodes n — 50 is 
significantly greater than the maximal number of indices k — 10. 
As usually log-uniform grid is used which has more nodes at low 
altitudes. 

The system is solved by S VD with regularization of the solu- 
tion by discarding small singular values s, < 10~ 3 So- The quality 
of the solution was controlled by two parameters: maximum de- 
viation from 1 (not accounting the initial part of the curve) and 



Table A.l. Coefficients cj (10~ 15 m 7 ^ 3 ) of h decomposition in 
WFs U(h) when we used all indices, 4 normal indices, 4 normal 
indices for original MASS device — columns 2, 3, 4 respectively 



Aperture 


c j 




c j 


A 


3.229 


2.981 


2.504 


B 


3.191 


-3.641 


-2.823 


C 


0.080 


2.880 


2.960 


D 


-0.826 


0.273 


-0.730 


AB 


-5.040 






AC 


1.124 






AD 


0.013 






BC 


-0.810 






BD 


0.547 






CD 


0.166 






Noise Factor 


4.16 


2.21 


2.54 



noise enhancement factor NF = Q^jCj) ' 2 / T*j c j- The solution 
depends on how the system was weighted. We used an implicit 
weight by grid density and some explicit weight. 

Without explicit weight, the solution has the first maximum 
at the minimum altitude but strongly oscillates and has large 
NF > 20. The study of additional weight in the form of ~ h p 
has shown that if p « 1 the solution is good enough: it is close 
to 1 and has moderate NF value. 

The approximating curves Au(h) calculated with all 10 in- 
dices and only 4 normal indices are shown in Fig. |4] The curve 
for aperture A for which an index and temporal covariances 
are used in the DESI method to assess the atmospheric coher- 
ence time is also shown. The coefficients cj are presented in 
Table I A. 1 1 for all variants. They were calculated for a typical 
set of MASS/DIMM apertures, its detector response and AO V 
spectral class of light source. 

The approximation which uses all indices is close to 1 for 
all nodes of the typical grid for restoration of OT profile. There 
is a slight (« 0.06) excess at the altitude of 0.7 km and the dif- 
ference for other nodes is lower. The approximation with only 4 
normal index looses 25% of the contribution of 1 km layer and 
almost 75% of the contribution of 0.5 km layer. It is expected 
that the decomposition in the complete set of indices will give 
an estimation of to with smaller systematic errors because the 
noise properties of cross-indices are not worse than that of nor- 
mal indices and the factor NF in the first variant is only twice as 
large. 
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